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agreement with experimental data is reasonable although it seems very difficult to obtain a perfect 
description in any case. 

\ PACS numbers: 12.39.Pn,13.40.Dk,14.20.-c,14.40.-n 

1 Keywords: Potential models; Electromagnetic mass differences; Baryons (including antiparticles); Mesons 

o 

(N 



I. INTRODUCTION 



Quantum Chromodynamics (QCD) is believed to be the good theory of strong interaction. It has met with numerous 
successes in many domains. However, in the low energy regime, it is extremely difficult to handle because of its non 
C^l perturbative character. Lattice calculations become more and more reliable but still remain very cumbersome, time 
K*" consuming and not always transparent for the underlying physics. This explains why, in the meson and baryon sectors, 
' a number of alternative simpler models were introduced. Among them, the non relativistic quark model (NRQM) is 
, very appealing because of its high simplicity, its ability to treat properly the center of mass motion, and the large 
• number of observables that can be described within its framework. 

\ In NRQM formalism, the dynamical equation is the usual Schrodinger equation including a non relativistic kinetic 
' energy term plus a potential term Q . There exist a lot of different numerical algorithms to solve the two-body and 
^> ' three-body problems with a good accuracy (see for instance Refs. @,@il3l)- 

(-H , Nowadays, it becomes more and more frequent to use a relativistic expression for the kinetic energy operator. The 
Q_i' resulting dynamical equation is known as a spinless Salpeter equation. It has several advantages as compared to the 
' Schrodinger equation 5, 6] and the corresponding numerical algorithms are now well under control (see for instance 

Refs. ^^JtiSiiM)- This kind of models and NRQM both use potentials, and are called potential models. 
(-H In potential models (and indeed in many other QCD inspired models) , the quark degrees of freedom are no longer 

• • the bare quarks of the QCD lagrangian, but are quasiparticles dressed by a gluon cloud and quark-antiquark virtual 
. !^ [ pairs. They are called constituent quarks and the most visible modification is the necessity to use in potential models 
I a quark mass substantially larger than the bare mass. In principle the bare quark-quark potential should also be 
}—{ ' modified and folded with the quark color density to give the final potential to be used in a Schrodinger or in a 
5^ , spinless Sal pet er equation. In recent quark-quark potentials appearing on the market, this effect is taken into account 
0,0,0,0]. Actually, this effect has been already considered since a rather long time (see i.e. Ref. 0,0|). The 
resulting spectra of these models are in correct agreement with experimental data. However, it seems very difficult to 
get, in a unified scheme (same form of the potential and same set of parameters), a g ood description of both meson 
and baryon properties [la.ll7j. For example in the well known works of Isgur [l4l Il5j |. the string tension is different 
in the mesonic and in the baryonic sector. Recently, progresses in obtaining a unified description have been achieved 
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The spectra are only a part of interesting observables, and the vaUdity of a model should be tested on other 
observables, especially if they are very sensitive to the form of the wave function. Electromagnetic properties are 
best suited for such a study, because the basic QED formalism is very well known and precise, and thus the possible 
uncertainties coming from mechanisms or wave function (itself depending on the much less known strong interactions) 
are more conveniently identified. 

In this paper, we focus our study on the electromagnetic splitting between charged hadrons, both in the meson 
and baryon sectors, within the framework of potential models. Since the earliest works on charmed mesons jl9l |. and 
baryons ia, a number of similar studies were performed in the past P, 0, |^ . Essentially three different 

sources for the splitting were identified: a small mass difference between up and down quarks, the coulomb interaction 
between charged quarks, and their dipole-dipole interaction p^ . All of them seem to have an important effect and 
the final result is a very subtle interplay among them. This explains why a very proper and precise treatment must be 
invoked, and also why this observable is very interesting. The current mass ratio for the u and d quarks is probably 
comprised between 0.2 and 0.8 ^2d\. However their absolute values are presumably weak (a few MeV) with respect to 
the QCD scale parameter A. Thus the spontaneous symmetry breaking induces large constituent masses for the u and 
d quarks; the corresponding values are very close (making the SU(2) isospin symmetry rather good) but nevertheless 
different. This small difference is a first source for the isospin splitting. Moreover, quarks being charged particles, 
the coulomb interaction is obviously present (very often, it has been treated as a perturbation). The dipole-dipole 
interaction (or hyperfine interaction) is a consequence of relativistic corrections to the coulomb potential. 

The same problem was also undertaken using formalisms relying more basically on fundamental QCD. In Ref. 
|25l |. an heavy quark effective theory (first order in l/mg) is adopted to study the splitting in heavy mesons, using 
dispersion relations. In Ref. a chiral field theory is employed to study several splittings in ordinary and strange 
sectors. In Ref. [l^, the authors are interested in the charmed and bottom meson sectors with a formalism based 
on Cottingham formula. In Ref. a tadpole term is introduced to deal with the u — d mass difference in some 
hadrons. But in most of these studies, the authors limit themselves to very restricted samples, either in meson or in 
baryon sector. We think that so few states to test the validity of a model can be questionable. 

In our paper we want to deal with all the known splittings both in mesons and in baryons in a consistent approach, 
and to push the potential model study further in several domains. First we want to perform a precise and complete 
treatment, avoiding perturbative expressions as much as possible. Second, we introduce the "contact term" , which 
arises on equal footing as the dipole-dipole relativistic correction, but which is neglected by most authors [l|, |23| ■ 

Finally, our most important improvement to our point of view is the use of a dressed electromagnetic interaction 
between quarks. Since the constituent quarks are quasiparticles, the electromagnetic interaction should also be 
modified as compared to the bare one, in a very similar way to the quark-quark strong potential. However the 
electromagnetic lagrangian is different from the QCD lagrangian and the electromagnetic density for the quark, playing 
a role in the splitting, has no reason to be identical to the color density occurring in the quark-quark strong potential. 
Such an approach has already been proposed in Refs. 0,|23| but with a different form for the electromagnetic density. 
Moreover, in these works, different sets of parameters have been used for mesonic and baryonic sectors separately. 

In order to see the sensitivity of the results on the treatment of the strong interactions for the quark dynamics, we 
investigate the splitting produced with two types of wave function, one resulting from a phcnomcnological non rela- 
tivistic hamiltonian (ALl) 29, 30] and another with a semi-relativistic hamiltonian (called here BSS) 18]. Moreover, 
since our aim is to consider mesons and baryons on equal footing, it is important to consider interquark potentials 
that lead to a correct description of both sectors. This is rather difficult to encounter. Both hamiltonians considered 
here are suited for that. The first potential (ALl) relies on the so-called funnel or Cornell potential [sil Is^. It is 
completely phenomenological and one can consider that the dressing of the quarks is included and simulated in the 
value of the various parameters. The second one (BSS) starts with more fundamental QCD grounds and is based on 
instanton induced effects 0,|13> dressing is explicitly taken into account but there remain nevertheless some free 
parameters that are adjusted on the spectra. 

The paper is organized as follows. The next section presents in more details the strong potentials and the way to 
solve the two-body and three-body problems. The third section deals with the electromagnetic interaction responsible 
of the splitting. The results of our calculations are presented and compared to data in the fourth section. Conclusions 
are drawn in the last section. 



II. POTENTIALS AND WAVE FUNCTION 
A. Strong potentials 

As stated in the introduction, we use in this paper two kinds of interquark potentials: one that must be used in a 
Schrodinger equation (ALl) and the other in a spinless Salpeter equation (BSS). They both depend on the relative 



distance r between the interacting quarks and are able to describe in a satisfactory way both the meson and the 
baryon spectra. However the BSS potential is suited, because of its underlying QCD basis, only for the light quark 
sectors (u, d, s quarks). They differ by the type of kinematics and by the manner to deal with spin splitting; although 
they give spectra of similar good quality, the corresponding wave functions can differ appreciably. Their derivation 
and their parameters have been reported elsewhere, and here we just want to point out the essential features and 
stress their differences more explicitly. In both models, the u and d quarks are assumed to have the same mass. In 
the following, they will be noted by the symbol n (normal or non strange). 



1. ALl Potential 

The ALl potential |3fl| |. developed for a non relativistic kinematics, contains the minimum ingredients necessary to 
get an overall reasonable description of hadronic resonances, a central part Vc and a hyperfinc term Vh 

V.,{r) - -^A, . A, [V^'^\r) + T/f \r)] . (1) 

The central part is merely the Cornell potential composed of a short range coulombic part, simulating the one-gluon 
exchange mechanism, and a long range linear term, responsible for the confinement (an additional constant is very 
important to get the good absolute values) 

V^''\r)^-^ + ar + C. (2) 

The color dependence through the Gell-Mann matrices A in relation comes from one gluon exchange. There is no 
reason, except simplicity, that such a structure is kept for the confining and constant parts. Nevertheless, this ansatz 
works well for both meson and baryon sectors. 

The hyperfine term has a short range behavior and is chosen as a gaussian function 

V^fe)^ - Stt exp(-rV4.) 

The interesting property, as compared to Bhaduri's [s^ or Cornell's 31] potentials, is that the range atj of that force 
does depend on the flavor 

, / 2m,m,- \ ^ 

a,, = A ^ . (4) 
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A kind of dressing is realized with the gaussian function, since the theory predicts a delta contribution for the hyperfine 
potential. The various parameters, including constituent masses, have been determined on the spectra by a best fit 
procedure. Although very simple, this potential does already a good job in hadronic spectroscopy. 

The ALl potential depends on quark masses in the strong hyperfine term. So a variation of the quark masses, as 
necessary for treating the electromagnetic splitting, has a strong effect on the meson masses. We have checked that 
this may induce wrong results in the case of a perturbative calculation. 



2. BSS Potential 

The BSS potential [l^ . developed for a semi relativistic kinematics, contains the Cornell potential and an instanton 
induced interaction. The Cornell potential has the same form as in the ALl potential, but the constant interaction is 
different for meson and baryon sectors. Contrary to the ALl potential, the instanton formalism is based on a SU(3) 
flavor symmetry and can hardly be generalized to the heavy quark sector. 

The instanton induced interaction provides a suitable formalism to reproduce well the spectrum of the pseudoscalar 
mesons (and to explain the masses of r/ and rj' mesons). The interaction between one quark and one antiquark in a 
meson is vanishing for L ^ or S ^ states. For L = 5 = 0, its form depends on the isospin of the qq pair 

• For / = 1: 

Vi{r)^-8gSir); (5) 



• For / = 1/2: 



Vi{r) = 



-8<?'5(f); 



(6) 



• 7 = 0: 



^.M-8(^^, ^^'M(f), (7) 



The para meters q and g' are two dimensioned coupling constants. Between two quarks in a baryon, this interaction 
is written 



in the flavor space (1/a/2(|uu) + \dd)), \ss)). 
mete rs q 

mill 

Vi (r) = -4 (.gP[""l + 5'p["'*l ) P^=°5{r) , (8) 

where P^='^ is the projector on spin 0, and P^'''' ' is the projector on antisymmetrical flavor state qq' . We have shown 
in Ref. |l8j that this model is able to describe correctly meson and baryon spectra in a consistent way if the following 
interaction is added for baryons only 

where Cj is a new constant. 

The quark masses used in this model are the constituent masses and not the current ones. It is then natural to 
suppose that a quark is not a pure pointlike particle, but an effective degree of freedom which is dressed by the 
gluon and quark-antiquark pair clouds. The form that we retain for the color charge density of a quark is a Gaussian 
function 

^W-(Tv^^-p(-'-Vr^). (10) 

It is generally assumed that the quark size T depends on the flavor. So, we consider two size parameters r„ and Tg 
for n and s quarks respectively. More details are given in Refs. hj,h3lhj. 

The instanton interaction acts differently on the symmetric {ud + du) and antisymmetric (ud — du) flavor states; in 
a framework where m„ — m^, the isospin formalism can be introduced to classify the basis states and this property 
can be taken into account without problem. As soon as m„ ^ nid, the up and down quark must be considered as 
different and the instanton becomes very difficult to handle. Thus, for baryons, the three-body treatment with BSS 
is done perturbatively. In this case, the potential has no mass dependence and the problem mentioned above for the 
ALl potential does not appear. 



B. Numerical Techniques 

There exist many numerical algorithms to compute the radial function of a meson. In this paper, we use the method 
based on Lagrange mesh, which is very simple, very precise and very fast, and for which a recent work has shown that 
relativistic kinematics can be handled without any problem [^. Thus the same method can be applied for both ALl 
and BSS potentials. To get a very high relative precision, around 10~^°, a typical number of basis states is iV = 60. 
Technical details can be found in Ref. 

Our method to solve the three-body problem is based on an expansion of the space wave function in terms of 
harmonic oscillator functions with two different sizes. This method, which is an old one (37j . was given up during a 
long time and renewed recently in Ref. Q (many authors use a numerical treatment based on harmonic oscillators 
but with a single variational parameter, for instance in Ref. [l5l|V In this paper, it was shown that the precision 
achieved was similar to the stochastic variational method Q. This last method was also shown to be competitive 
with more conventional methods, such as Faddeev formalism |38| |. The number of quanta is the relevant quantity for 
convergence; a very good relative precision of the order of 10~^, necessary in our study, is achieved if we include all 
basis states up to 20 quanta. More details on the method can be found in Ref. 0- 



III. ELECTROMAGNETIC INTERACTION 



A. Bare potential 

The electromagnetic interaction between two pointlike particles i and j of charges Qi, Qj and masses rrii, rrij is 
very well known. In addition to the usual coulomb potential ?7coui, relativistic corrections (at lowest order) give rise 
to contact, hyperfine, tensor, symmetric and antisymmetric spin-orbit, and Darwin terms. Tensor, spin-orbit, and 
Darwin terms are complicated and presumably their effects are weak, so that people neglect them. The hyperfine 
interaction Uhyp does play an important role and is included in all serious calculations. The contact term ?7cont 
is usually discarded with no other justification than simplicity. In this study we keep it, so that our total bare 
electromagnetic potential writes 

C/W(r) ^ ([/_,)!? (r) + (t/co„t)!?(r) + {U^yJ^\r), (11) 

with 



(?7coui)!-^(r) = Q^Q,-, (12) 
(C/cont)iy (r) = -|g.Q, (^-^ + aSir), (13) 

(U.y^^(r) = -|^«<5(r)s..s„ (14) 

where a is the fine structure constant. Let us stress a point that deserves further discussion. One can wonder why 
such a contact term is not kept as well in the strong potential ALl. This term comes from one gluon exchange and is 
proportional to as] but this parameter is essentially phenomenological and the Cornell form of the central potential 
is a kind of averaging of all unknown effects coming from QCD, including the contact term. Our philosophy is to 
maintain as far possible the form of the strong potential, that was shown to do a good job in hadronic physics. On 
the contrary, the situation is completely different in the electromagnetic hamiltonian Hll|l because the corresponding 
expression is the exact one and does not suffer from ambiguities or phenomenological uncertainties. This is why we 
decided to include the contact term, which is usually neglected, in order to grasp quantitatively its effect and to justify 
a posteriori the vahdity (or not) of neglecting it. 

A perturbative treatment is compulsory for such a bare electromagnetic potential, owing to the presence of the 
delta function. It was shown that, in some cases, this a very bad approximation [12j |. 



B. Electromagnetic quark density 

One way to introduce the very complicated mechanisms that transform the bare quarks (pointlike) to constituent 
quarks is to assume the existence of a phenomenological density distribution for a quark. The electromagnetic density 
should in principle depend on the internal structure of the system and contain a lot of unknown processes (including 
probably interferences between QCD and QED) leading to isospin symmetry breaking. 

This means that a constituent quark located at r is generated by a bare quark located in r' with a certain probability 
distribution (or density) p(r — r'). To have an appealing physical meaning this density must be a peaked function 
around zero with a certain size parameter, that depends in principle on the quark flavor. Moreover we require the 
natural property that, for a vanishing size, the constituent quark becomes pointlike and is identified to the bare quark. 
Mathematically, this means that the limit of the density p(u) for a vanishing size is the delta function S{u). Another 
natural property is that the density is isotropic. The last required property is that the integral of the density over 
the whole space is unity. 

The most popular densities are of lorentzian, gaussian or Yukawa type. For instance, the authors of Refs. pll |23 | 
adopted a gaussian form. In this study we choose an electromagnetic Yukawa density. There is a precise reason 
for that: this density is the leading ingredient of the meson charge form factor. It is an experimental fact that the 
data accommodate rather nicely to a Yukawa density, giving a form factor with an asymptotic behavior [sMIs^. 
instead of a gaussian density which leads to gaussian asymptotic behavior. Keeping in mind the previous remarks, 
the adopted density for quark of flavor i is 



Pii^) = 1 2 ' 



(15) 



where 7^ is the electromagnetic size parameter. 

The dressed potential U is obtained from the bare potential [/^''•' by a double convolution over the densities of each 
interacting quark 

f/„- (r) = Jdudw (u) p, (v) U^f (r + V - u) . (16) 
With a trivial change of variable, it is quite easy to transform this double folding into a single one 

f/.,(r)= j dr' (r')p,:,(r-r'), (17) 

with the definition of the new density 

Ptj{'^) = J dvPtMPji'^-'^)- (18) 
Applying Eq. (|18|l to the Yukawa density 115|l , one gets 

P^A^)^^e--^" if 7.-7, (19) 

and 



'«(")-3;hF^(— -— ) <») 

= (^2 ^^2) (7fp»(u)-7>j(u)) if J^^lr (21) 
C. Dressed potential 

Starting from Eq. Hll|l and using Eqs. H17() and H19(l . it is easy to obtain analytically the dressed potential 

Uijir) = ([/coui)y(r-) + {Uconthir) + (C/hyp)y(r), (22) 
in the case of two interacting constituent quarks with the same size 7^ = jj — 7. It writes 



iUcoui)ijir) = aQiQj 



1 / / \ e^'"/''' 



(23) 



(f/„,,),W^-^|ii(ij + ij).-A (24) 

Doing the same thing with Eq. (|21|l , one obtains the electromagnetic dressed potential in the case of different inter- 
acting quarks 

(t/couO.. M = - - ^2!: ^2 , + ^2^2 , - (26) 



aQ^QJ fj_^ J_\ e-'-/T' - e-'-/T^ 
8(7,— 7j^) \m^i 



(U.oMr) = -^7^r%T + ^ (27) 



(^hyp).,(r) = 2, - ^ (28) 

3mimj(7f-7|) r 

It is possible to treat the electromagnetic potential [/ as a perturbation. This procedure is only used for the baryons 
with the BSS potential (see Sec. Ill A 2)) . 



IV. RESULTS 



A. Determination of the parameters 

We do not want to introduce a lot of new parameters; here we restrict the number of free parameters to the minimum 
unavoidable. In particular the parameters of the strong potential are maintained without modification. Important 
information is contained in the electromagnetic size of the quarks 7^. This size could depend on the flavor and on the 
electrical charge of the quark. If the dependence on charge were dominant, we could expect that sizes of quarks with 
the same charge are similar. Preliminary calculations have shown that the size must be strongly reduced for heavy 
quarks, showing that the dependence on mass must be dominant. Consequently, in order to restrict again the number 
of parameters, we assume that the sizes of the u and d quarks are the same. Thus we impose 7u = 7ri = 7n- For ALl 
and BSS models, we have 7„ and 7s as free parameters. An observable that is very sensitive to those parameters is 
the charge mean square radius. 

In the mesonic sector, the bare charge square radius operator for pointlike quarks is defined by 

(r^)('')=^e.(r.-R)^ (29) 

where is the charge for quark i, its position, and R the position of the center of mass. For constituent quarks, 
this expression has to be folded with quark density and should be written instead 

2 

(r^) = E - ^^('•' - (30) 

i—l 

Averaging quantity (|30|l on the meson wave function provides us with the charge mean square radius of the meson. 

Performing the calculation, one finds that this observable is a sum of a term, that can be called the bare radius <^r^) 
(which is essentially the mean value of quantity H29|l on meson wave function), plus a term which is essentially the 
dipole moment of the density. With a Yukawa density, one has explicitly 

(r^)=(r2)^n6^e.7f. (31) 

i=l 

The dynamical contribution to the square radius is entirely contained in the bare quantity, whose expression is very 
well known and is not recalled here. 

From relation (|31|l . one sees that the pion radius depends only on 7„ and it is used to determine this parameter. 
The kaon radius depends on 7„ and 7s; since 7„ has already be determined from the pion, the kaon radius is used to 
determine 7s. In addition, the ALl potential needs the further determination of 7c and 7b. Since the radii for D and B 
resonances are not known experimentally, we determine the corresponding sizes by requiring a smooth behavior versus 
the mass. Moreover, the uncertainty on the kaon radius is rather large. Fortunately, we checked that the isospin 
splitting does not depend too much on the precise value of the quark size. Our accepted values for the electromagnetic 
sizes are summarized on tabled] 

Other parameters that need to be changed are the quark masses. In fact, the only important ingredient for the 
splittings is the mass difference between the down and up quarks: A = — m„. In view of this, we choose to 
maintain the s, c, and b quark masses at their non perturbative value fhi = rrii, and to keep the average value of 
the isospin doublet at its non perturbative value (m^ + rh^jl = rrin- The size parameters being determined once 
and for all on charge radii, we have only one free parameter A at our disposal to try to reproduce all the known 
electromagnetic splittings. To really see that this is a very big constraint, let us recall that for doing the same job, 
authors of Refs. |2^|23| used 4 free parameters and Genovese et al 2 free parameters. The 7r° and p° resonances are 
considered here as nn systems, composed of fictitious quark and antiquark of mass m„ = (to„ + md)/2 = m„. Doing 
this, we neglect mixing of J = and I ~ 1 components; this is a good approximation up to second order in A/m„. 
One can imagine several strategies to determine the parameter A. We first remarked that if we fit A on the mesons, 
the baryons were very badly reproduced, while if we fit it on the baryons, the mesons were spoiled less dramatically. 
Moreover, among the baryons, some splittings are more affected than others by a small change of A. Thus, we decided 
to fit this parameter on one of the most sensitive and well known splitting, namely — 'S'^ — 4.807 ± 0.04 MeV. 

One aim of our study is to see the respective influence of each component of the electromagnetic potential H22|l . 
The coulomb part seems to us unavoidable, so we will consider four different approximations in the following: 

• C: Electromagnetic potential restricted to the coulomb term (I23f) or H2t)|l alone; 



• CC: Electromagnetic potential restricted to coulomb and contact H24() or (|27|l terms alone; 



• CH: Electromagnetic potential restricted to coulomb and hyperfine H25|) or ()28|l terms alone; 

• T: Total electromagnetic potential H22() taken into account. 

Each approximation requires its own A parameter, but the 7^ parameters can be maintained to their values of TableQ] 
The corresponding results are gathered in Table ^ It is worth noting that the study of any of these approximations 
is a complete calculation by its own and needs a computational effort identical to the total treatment. 

It is funny that for both potentials, increasing values of A are obtained with approximations CH, C, T, CC 
respectively. Is it a property independent of the strong potential (realistic enough to reproduce baryon and meson 
spectra)? We have no answer. 

B. Experimental data 

A number of experimental data, concerning the electromagnetic splittings, exist for both the mesonic and baryonic 
sectors, and for both light quarks (u, d, and s) systems and systems containing at least a heavy quark (c and b). 
The use of ALl potential allows to study all data, while BSS is restricted to the light quark domain for the reasons 
mentioned in Sec. Ill A 21 The values for the splittings are of order of some MeV and some of them are known with 
good accuracy. 

From the splitting in the nucleon case and many others, the d quark mass is presumably larger than the u quark 
mass. From a naive argument based on rest masses only, the hierarchy of the splittings can be explained although the 
quantitative value needs much refined explanation. However, contrary to this naive argument, a number of puzzling 
questions still arise and some of them are not solved in a satisfactory manner by the up to date theoretical studies. 
Let us list some of them: 

9 n — p — 1.293 MeV is much weaker than 7r+ — 7r° = 4.594 MeV, despite the fact that a very naive quark model 
gives n{udd) — p{duu) = A and ^^{ud) — Tr^{[uu — dd]/\/2) = 0; 

• 7r+ — TT*^ has a large positive value while — has a small negative value (a positive value is compatible with 
the error bar but with a very low value in any case); 

• D2 — D2 ~ (or even may be negative) while the values for other cfi states are largely positive; 

• S^""*" is the highest level of the multiplet, in contradiction to what one expects from naive arguments where it 
should be the lowest one; 

• All the members of the charmed 5 verify > S+ except the particular states S|^* < S+* which satisfy the 
opposite relation. 

These remarks illustrate the fact that the mass difference between down and up quarks cannot be the only - or even 
leading - ingredient to explain the splitting; the internal dynamics also plays an important role. An important aim 
of this paper is devoted to examine this question. 

C. Influence of various approximations 

In this part we want to study the effect of using various approximations C, CC, CH and T of the electromagnetic 
hamiltonian. An exact treatment is performed, except for the baryons with BSS potential, as explained previously. 

Let us first present the results obtained with ALl potential. They are gathered on Table UTTl for the mesons and on 
Table HVl for the baryons. Few comments are in order: 

• The various approximations differ significantly one from the other. This proves again that the electromagnetic 
splitting is an observable very sensitive to the physical content put in it. A more quantitative comparison is 
relegated later on. 

• Although the experimental data cannot be reproduced with a good precision, the calculated values have the 
good order of magnitude and respect more or less the hierarchy. We mean by this that the order of a given 
multiplet is generally the good one and that a large (small) theoretical splitting corresponds to a large (small) 
experimental one. Let us recall that we have only one free parameter A, which has been fitted on the — E*^ 
value. 



• From time to time the sign is wrong (order is opposite to the experimental one) , but this effect occurs generaUy 
when compatibihty is not excluded due to error bars, or at least when the experimental uncertainty is large. 

• Let us comment on the puzzling questions raised in the introduction of this section. The splitting among the 
pions is larger than the splitting among the nucleons, but not enough to claim that the problem is solved. The 
/)+ — mass difference is still wrong but the value is lower than the naive expected one. The characteristic 
pattern for the D mesons is good but the quantitative values fail. The order in the multiplet is a good one 
and the quantitative value are also very good; this case is really a success (let us mention that in Ref. ,22j the 
theoretical value agrees well with the experimental data, but since then the data changed and the agreement 
is less good with the new value). The problem for the Sc is not solved but the order is the good one for 3 
resonances over 4. 

Let us now have a look on the situation concerning the wave functions arising from the BSS potential. The mesonic 
sector is presented in Table fvl and the baryonic sector in Table IVll The same comments as for ALl can be made for 
BSS generally, but in this case the puzzle concerning the p has found a solution. Unfortunately, the nucleon splitting 
is much too large. 

Moreover the BSS results are quite different from ALl results, indicating that, in order to describe the electromag- 
netic splitting, not only a good form for the electromagnetic part of the interaction is important, but also the strong 
one via the wave function. This conclusion is very important. In order to draw this conclusion, it was necessary to 
use very different strong potentials that lead equally good results on spectra both in the mesonic and in the baryonic 
sectors. These requirements are not easy to be met. ALl and BSS allow such a fruitful comparison. 

In order to grasp more quantitatively the influence of each approximation, as well as to compare the effect of 
the strong potential, we calculate a chi-square value on the experimental sample (in fact a chi-square divided by 
the number of data in the set). In order to avoid a too much important weight on very precise values, a minimum 
uncertainty at 0.1 MeV has been assigned arbitrarily to those values. The results are gathered in Table IVTll Let us 
note that this chi-square has no precise statistical relevance; it is just a convenient way to get a synthetic view of the 
comparison between the various approximations. To have a more refined analysis, we separate the sample for meson 
(denoted by "M" in the Table), for baryon ("B" in the Table), or the entire sample of hadrons ("H" = "M"-h"B" in 
the Table). Moreover we also distinguish between light sector (u, c?, s denoted by "1" in the Table), heavy sector (c, 
h denoted by "h" in the Table), or all sectors (denoted by "a" = "l"-t-"h" in the Table). For instance the line "hB" 
means a chi-square calculated on heavy baryons, "IM" calculated on light mesons, "aH" on the entire sample, . . . 

Interesting remarks can be emphasized: 

• The baryonic sector is explained in a much more satisfactory way, whatever the chosen strong potential. This 
is the consequence of our arbitrary choice to adjust the free parameter A on a baryon resonance. Should we 
have chosen to fit A on a meson resonance, the mesonic sector would have been much best reproduced, but at 
the price of a dramatic spoiling of the baryonic sector, and with an overall worse description. It appears that a 
consistent description of electromagnetic splittings for both mesons and baryons is hardly feasible with a model 
containing only one free parameter. Let us stress that our model probably do not contain all possible sources of 
splitting. For instance, the pion mass difference can be explained with the vector meson dominance assumption 
|39j. Nevertheless, this last model relies on pointlike pion while mesons are composite particles in our approach. 
So it is difficult to compare such so different phenomenological models of QCD. 

• Concerning the BSS calculation, the exact formalism (T approximation) is the best one for mesons, but it is the 
CH approximation that is the best for baryons, and the C approximation (presumably the crudest one) for the 
entire set. Anyhow, addition of the contact term deteriorates seriously the results. 

• Concerning the ALl calculation, several conclusions are drawn. For baryons, all approximations are roughly of 
the same quality. What is gained on data by one approximation versus the others is lost on another data. But 
curiously CC is the best approximations and T the worse one, so that the conclusion is rather different from 
BSS case. For mesons, CH is always the best approximation and CC the worse. For the totality of the sample 
this conclusion remains, while the exact treatment (T) is just a bit better than the crudest one (C). Here again 
the contact term has a very bad influence. 

• ALl and BSS calculations give results of comparable quality for the whole set, but BSS is much better for mesons, 
and ALl much better for baryons. This may be explained by the fact that only a perturbative treatment (and 
not an exact one as in ALl) can be performed with BSS in the baryonic sector. 

• It is interesting to compare our results with some previous studies treating both meson and baryon electromag- 
netic splittings: 



— In Refs. |2lj,|22|, the results look at least as good as ours, but we want address two comments. First, some 
data have changed since that time (for instance, the old value of 1.8 MeV for the Sji"^ — S° splitting is 
now 0.35 MeV) and new ones are now available (for instance, their calculated value of 4.4 MeV for the 
D2 — D2 is now experimentally estimated at 0.1 MeV). Second, the value of the electromagnetic quark 
size is determined by a best fit on the splittings, while in our case, the same parameter is fixed by meson 
form factors and cannot be considered as a free parameter. 

— In Rcf. 23], the results look also rather good, but there exist several differences with the present work. 
First, the contact term is absent. Second, only the hyperfine term is dressed in order to avoid a collapse. 
Third, two free parameters were used, the u and d masses separately. 

No doubt that if we have allowed more free parameters in our model (for example releasing the constraint (m„ + 
fhd)/2 — vrin), our results would have appeared better. But this was not the ultimate goal of our paper; we were 
interested in discovering what are the necessary ingredients to explain the splittings. 

Another important point for the consistency of our approach is that, with the new values of the u and d masses, 
the absolute masses for the meson and baryon resonances are also well reproduced. This point is not often studied 
in previous works. Just to convince the reader that our formalism is able to provide a good description of meson 
and baryons simultaneously, we present below the absolute masses of one member of the niultiplet (the other can be 
obtained using the values of the splittings given above) for both mesons (Table EnH) and baryons (Table |Tx|). The T 
approximation is used and the treatment is exact except for baryons with the BSS potential. 

The small discrepancy for light baryons in the case of ALl potential can be attributed to a three-body force [s^ . 
which is mass dependent. The instanton does not give rise to a three-body force for the baryon 0, and the 
parameters have been fitted directly to the absolute masses of mesons and baryons. The agreement with experimental 
data is rather satisfactory. 



V. CONCLUSIONS 



In this paper, we calculated the electromagnetic splitting on hadronic systems. This observable is a very sensitive 
test of the formalism because it results from a very subtle and fine balance between several physical ingredients. In 
order to concentrate on the physical aspects of the problem, we were very cautious in the numerical treatment, both 
for the two-body and three-body problems. Thus we are very confident with our numerical results, and interesting 
conclusions can be emphasized. 

As compared to previous works we considered the total electromagnetic hamiltonian (excepted Darwin, spin-orbit, 
and tensor forces that we believe to play a very minor role) . In particular we took into account the so-called contact 
term that is usually neglected. Moreover, we treated the quarks as constituent particles with an electromagnetic 
size that modifies the form of the electromagnetic interaction. Lastly, we based our calculations on wave functions 
resulting from two different strong potentials: the phenomenological ALl potential to be used with a non relativistic 
kinematics energy operator and the more fundamental BSS potential, including instanton effects, to be used with a 
relativistic kinetic energy operator. The size of the quarks were determined in order to reproduce the pion and kaon 
charge form factors. Only one free parameter, the mass difference between down and up quarks, is left at our disposal; 
it was fitted on the — TP splitting. We want to emphasize that in this paper, with only one parameter, we try to 
reproduce the totality of the experimental data, including mesons and baryons. To our opinion, this is a condition 
to pretend to some consistency in the formalism. This condition is rarely met in previous calculations, since authors 
very often restricted themselves either to mesons (or even to less restrictive domains) or to baryons. 

We first showed that dressing the strong and electromagnetic interaction is a necessity to obtain in a consistent 
way both the hadron spectra, the meson form factors and the electromagnetic splitting. 

By comparison of ALl and BSS results, we stressed that the strong potential, via the wave function, is an important 
ingredient in the description of electromagnetic splittings. In our particular case, ALl gives a better description of 
baryons, and BSS a better description of mesons. 

But we proved also that the electromagnetic hamiltonian is equally important for explaining the splittings. Adding 
or removing a term (contact or hyperfine) has a non negligible influence. In particular taking into account the contact 
term spoils a lot the results, and curiously it is the approximation based on coulomb-|-hyperfinc (the usual ingredient 
for many people) which is globally the best. 

In our formalism the splittings are described in a reasonable way, specially owing to the fact that we have only one 
free parameter to try to reproduce 26 experimental data; the order among multiplet masses is generally correct and 
the values have the right order of magnitude. But the agreement is far from being perfect; some suggested puzzles 
have been solved but some others are still opened questions. New improvements must be done in future studies. 
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In 


7s 


7c 


7b 


ALl 
BSS 


1.225 
1.330 


0.200 
0.450 


0.040 


0.013 



TABLE L Electromagnetic sizes 7/ of the constituents quarks (in GeV for different flavors /, and for the two different 
strong potentials. 





C 


CC 


CH 


T 


ALl 
BSS 


13.2 
7.2 


14.4 
9.4 


12.6 
6.2 


13.6 
8.4 



TABLE IL Values of the parameter A = md — mu (in MeV), for the two potentials, and for the four different approximations 
presented in the text. 



Splitting 


Exp 


C 


CC 


CH 


T 


7T+ - vr" 


4.594 ± 0.001 


1.239 


0.828 


2.111 


1.693 




-0.5 ±0.7 


1.009 


0.857 


0.870 


0.713 


K°~K+ 


3.995 ± 0.034 


8.096 


9.290 


7.054 


8.109 


KO* - K+* 


6.7 ± 1.2 


1.132 


1.458 


1.139 


1.436 




6.8 ±2.8 


-0.794 


-0.793 


-0.763 


-0.757 


D+-D° 


4.78 ±0.10 


2.478 


1.977 


2.821 


2.308 


D+' - DO* 


2.6 ± 1.8 


1.513 


1.120 


1.428 


1.037 


Dt-D° 


6.8 ±5 


-2.038 


-2.409 


-1.839 


-2.161 


Dt-Dl 


0.1 ±4 


-2.058 


-2.411 


-1.937 


-2.235 


B° - B+ 


0.33 ±0.28 


-1.877 


-1.715 


-1.891 


-1.713 



TABLE III: Meson electromagnetic splittings (in MeV) calculated for 4 different approximations C, CC, CH, and T, as explained 
in the text. The wave functions result from the strong potential ALl and the numerical treatment is exact. The experimental 
values (Exp) come from Ref. [2^ . 



Splitting 


Exp 


C 


CC 


CH 


T 


n — p 


1.293 


0.89 


1.15 


0.91 


1.15 


A" - A++ 


2.25 ± 0.68 


2.68 


3.68 


2.79 


3.72 


A+ - A++ 


1.2 ±0.6 


0.57 


1.22 


0.73 


1.35 


S- -E" 


4.807 ± 0.04 


4.82 


4.82 


4.81 


4.76 


E" - E+ 


8.08 ± 0.08 


7.87 


8.41 


8.27 


8.55 


E-*-E0* 


2.0 ± 2.4 


3.26 


3.23 


3.04 


2.94 


E"* - E+* 


0±4 


1.71 


1.99 


1.68 


1.96 




6.48 ± 0.24 


7.12 


7.21 


7.38 


7.38 




3.2 ±0.6 


3.01 


2.91 


2.80 


2.66 


E^+ - EO 


0.35 ±0.18 


1.00 


-0.02 


1.35 


0.37 


E° - E^ 


0.9 ±0.4 


0.32 


0.64 


0.02 


0.33 




1.9 ±1.7 


1.37 


0.27 


1.33 


0.19 


=■0 = + 

' — 'C ' — 'C 


5.5 ±1.8 


2.81 


3.28 


3.01 


3.42 


=■0/ 

' — 'c ' — 'c 


~ 4.2 ± 3.5 


0.20 


0.49 


-0.08 


0.20 


' — 'c " ' — 'c 


~ 2.9 ±2.0 


-0.08 


-0.31 


-0.03 


-0.25 


' — 'C ' — 'C 


~ 4.1 ± 2.5 


3.09 


3.42 


3.24 


3.51 



TABLE IV: Same as Table 11111 for baryons. The theoretical uncertainty may affect only the last digit. 



Splitting 


Exp 


C 


CC 


CH 


T 




4.594 ± 0.001 


1.22 


-0.21 


4.11 


2.97 


p+-pO 


-0.5 ±0.7 


1.00 


0.40 


0.59 


-0.01 


K'>-K+ 


3.995 ± 0.034 


1.54 


3.68 


-0.97 


1.17 


K'"- ~ K+' 


6.7 ± 1.2 


2.88 


4.50 


2.63 


4.25 


K', - Kt 


6.8 ±2.8 


2.41 


3.49 


2.07 


2.83 



TABLE V: Meson electromagnetic splittings (in MeV) calculated for 4 different approximations C, CC, CH, and T, as explained 
in the text. The wave functions result from the strong potential BSS and the numerical treatment is exact. The experimental 
values (Exp) come from Ref. [2^ . 



Splitting 


Exp 


C 


CC 


CH 


T 


n — p 


1.293 


2.78 


4.28 


2.62 


4.08 


A° - A++ 


2.25 ± 0.68 


4.24 


8.21 


4.37 


8.34 


A+ - A++ 


1.2 ±0.6 


1.18 


3.69 


1.59 


4.09 


E- -EO 


4.807 ± 0.035 


4.83 


4.81 


4.81 


4.80 


E" -E+ 


8.08 ± 0.08 


7.68 


8.90 


8.47 


9.69 


E-* - EO* 


2.0 ±2.4 


4.98 


5.37 


4.06 


4.43 


E"* - E+* 


0±4 


3.12 


4.54 


2.88 


4.33 


^- .^0 


6.48 ± 0.24 


4.82 


4.50 


5.87 


5.55 




3.2 ±0.6 


5.09 


5.38 


4.17 


4.46 



TABLE VL Same as Table fvl for baryons. The theoretical uncertainty may affect only the last digit. Let us recall that, for 
technical reasons, the theoretical values were obtained with a perturbative treatment. 





ALl 


BSS 




C 


CC 


CH 


T 


C 


CC 


CH 


T 


IM 


568 


851 


317 


513 


350 


464 


500 


213 


IB 


3.3 


3.0 


3.8 


5.0 


34 


127 


23 


129 


IH 


205 


306 


116 


186 


147 


247 


194 


159 


hM 


119 


169 


90 


134 










hB 


3.0 


1.5 


5.9 


1.2 










hH 


51 


73 


41 


56 










aM 


344 


510 


204 


324 










aB 


3.2 


2.4 


4.7 


3.3 










aH 


134 


198 


81 


127 











TABLE VH: Chi-square values divided by the number of data in the sample. The meaning of approximations C, CC, CH, and 
T have been explained in the text. The meaning of the rows is related to the sub-samples taken into account: "1", "h", "a" for 
light, heavy and all sectors respectively; "M", "B", "H" for meson, baryon, hadron sectors respectively. 



System 


Exp 


ALl 


BSS 




134.98 


137.3 


145.8 


P" 


769.0 


769.7 


756.2 


K+ 


493.68 


486.6 


491.2 


K+* 


891.66 


902.7 


888.7 


Kt 


1425.6 


1333.4 


1377.5 


D+ 


1869.3 


1863.2 




D+' 


2010.0 


2016.4 




Dt 


2427 


2419.6 




Dt 


2459 


2451 




B+ 


5279 


5295 





TABLE VIII: Absolute masses for some mesons. Calculations are done with strong potentials ALl and BSS, and with the total 
electromagnetic hamiltonian. The experimental values (Exp) come from Ref. |2^ . 



System 


Exp 


ALl 


BSS 


V 


938 


994 


935 


A" 


1234 


1308 


1260 


A 


1116 


1149 


1105 


E" 


1197 


1233 


1201 


E-* 


1387 


1439 


1395 


H" 


1321 


1343 


1323 


H~* 


1535 


1560 


1522 




1672 


1675 


1646 




2285 


2290 


_ 




2453 


2466 


— 


^+ 

' — 'C 


2466 


2467 




^+/ 

' — 'C 


2574 


2572 


_ 


'— 'c 


2647 


2650 


_ 




2815 


2788 






2697 


2675 







5624 


5635 





TABLE IX: Same 



as Table IVTTTl for 



baryons. 



